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^ I Abstract 

"^ ! We present a traffic model that extends the hnear car-following model as well as the min- 

OC I plus traffic model (a model based on the min-plus algebra). A discrete-time car-dynamics 

describing the traffic on a 1-lane road without passing is interpreted as a dynamic program- 
ming equation of a stochastic optimal control problem of a Markov chain. This variational 
formulation permits to characterize the stability of the car-dynamics and to calculte the 
stationary regimes when they exist. The model is based on a piecewise linear approximation 
of the fundamental traffic diagram. 
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1. Introduction 
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\Q ', Car-following models are microscopic traffic models that describe the car-dynamics with 

stimulus-response equations expressing the drivers' behavior. Each driver responds, by choos- 
ing its speed or acceleration, to a given stimulus that can be composed of many factors such 
as inter-vehicular distances, relative velocities, instantaneous velocities, etc. We present in 
this article a car-following model that extends the linear car-following model p, l2l|, |lj, |l3 



as well as the min-plus traffic model [27 



The vehicular traffic on a 1-lane road without passing is described by discrete-time dy- 
/\ • namics, which are interpreted as dynamic programming equations associated to stochastic 

c^ ■ optimal control problems of Markov chains. The discrete-time variational formulation we 

make here is similar to the time-continuous one used by Daganzo and Geroliminis J7[ to show 
the existence of a concave macroscopic fundamental diagram. 

We are concerned in this article by microscopic traffic modeling with a Lagrangian de- 
scription of the traffic dynamics, where the function x{n,t), giving the position of car n at 
time t (or the cumulated distance traveled by a car n up to time t), is used. In the macro- 
scopic kinematic traffic modeling, Eulerian descriptions of the traffic dynamics are usually 
used, with the function n(t, x) giving the cumulated number of cars passing through position 
X up to time t (which coincides with the Moscowitz function |6| in the case of traffic without 
passing). The combination of a conservation law with an equilibrium law gives the well 



known ffist order traffic model of Lighthill, Whitham and Richards [26|, |29 . 
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In this introduction, we first notice that the same approach used in macroscopic traffic 
modehng, combining a conservation law with an equihbrium law, can be used to derive 
microscopic traffic models. By this, we introduce car-following models and give a review on 
basic and most important existing ones. In particular, we recall the linear and the min-plus 
car-following models, which are particular cases of the model we present here. Finally, we 
recall a theoretical result on non-expansive and connected dynamic systems which we need 
in our developments, and give the outline of this article. 

The ffist order partial derivative of the function x{n, t) in time, denoted x{n, t) expresses 
the velocity v{n,t) of car n at time t. The ffist order differentiation of x{n,t) in the car 
numbers Ax{n, t) expresses the inverse of the inter-vehicular distance y{n, t) = x{n — 1, t) — 
x{n,t) between cars n and n — 1 at time t. The equality of the second derivatives Ax{n,t) 
and Ax{n,t) gives then the following conservation law of distance. 

y{n,t) = -Av{n,t). (1) 

If we assume that a fundamental diagram Ve, giving the velocity f as a function of the 
inter- vehicular distance y {v = Ve{y)) at the stationary traffic, exists, and that the diagram 
Ve holds also on the transient traffic, then we have 

y{n,t)=v{n,t)/V:{y), (2) 

where V^{y) denotes the derivative of Ve with respect to y. 
By combining ([T]) and ([2]) we obtain the model: 

vin,t) = -VMn,t))Mn,t). (3) 

(|3]) is a car-following model that gives the acceleration of car n at time t as a response to 
a stimulus composed of the relative speed Av{n,t) and the term V^{y{n,t)). For example, if 
Ve{y) = voexp —a/y, where vq denotes the free (or desired) velocity, and a is a parameter, 
then V'(y) = aV{y)/y^ and (j3]) gives a particular case of the Gazis, Herman, and Rothery 
model [l3[. 

The simplest car-following model is the linear one, where the car dynamics are written 

Xn{t + T)=a{Xn-l{t)-Xn{t))+b, (4) 

where T is the reaction time, a is a sensitivity parameter, and 6 is a constant. That model 
can be derived from (^ with a linear fundamental diagram Ve- The stability of the linear 



car-following model 01]) and the existence of a stationary regime have been treated in [21 



Almost all car following models are based on the assumption of the existence of a be 



havioral law Ve- The latter has been taken linear in [17|, |2l|], logarithmic in [16|, exponential 



m 



28|, and with more complicated forms in other works. Bando et al. [3j used the sigmiudal 
function 

K(y) = tanh{y — h) + tanh{h), (5) 

where tank denotes the hyperbolic tangent function, and /i is a constant. 



Kerner and Konhauser[24|. Hermann and Kerner [22|, and then Lenz et al. [25|, and 



Hoogendoorn et al. [23| have used the function 

1,000 



Ve{y) = Vo 



1 + exp 



1-y 
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5.34 • 10' 



(6) 



where free velocity f n and 7 are parameters estimated from data. In [25|, 7 is taken equal 



to 7.5. See also 19, 20, 30 



The min-plus traffic model [27| is a microscopic discrete-time car-following model based 
on the min-plus algebra [2|. It consists in the following dynamics. 

(7) 



Xnit + 1) = min{x„(t) + vq, x„_i(t) - a}, 



where vq is the free velocity and a is a safety distance. 

The idea of the model ([7]) is that the car dynamics is linear in the min-plus algebra [2| , 
where the addition is the operation "min" and the product is the standard addition "+". 
In 271, the dynamic system ([7]) is written in min-plus notations x{t) = A® x{t — 1), where 
A is a min-plus matrix and ® is the min-plus product of matrices. It is then proved that the 
average growth rate vector per time unit of the system, defined by x = linii_j.oo x{t)/t satisfies 
X = iv,v,- ■ ■ ,v), where v is the unique min-plus eigenvalue of the matrix A, interpreted as 
the stationary car-speed. The fundamental diagram is then obtained 



V = mm(t>o,y — a) 
q = mm{vop,l - ap), 



(9) 



where y, q and p denote respectively the equilibrium inter-vehicular distance, car-flow and 
car-density. 

In ([7]), for t > 0, if x„_i(t) — x„(t) > vq + a, then the dynamics is x.„(t + 1) — x„(t) = vq. 
If x„_i(t) — Xn(t) < I'o + <^ then the dynamics is x„(t -|- 1) — x„(t) = x„_i(t) — x„(t) — a. 
Therefore, ([7]) is linear in both phases of free and congested traffic. The min-plus model ([7]) 
permits to distinguish two phases in which the traffic dynamics is linear, but with a sensitivity 
parameter (a in (jl])) equals to 1 for each phase. The linear model (jl]) is not constrained in 
the sensitivity parameter value, but it permits the modeling of only one traffic phase. The 
model we present here extends (jl]) and ([7]) in a way that an arbitrary number of traffic phases 
can be modeled, with flexibility in the sensitivity parameter value on each phase. 

For our model, we apply a similar but more general approach than the min-plus one used 
to analyze the dynamic system ([7]). Indeed, the dynamics ([7]) is additive homogeneous of 
degree oneu and is monotone q It is then non expansive U |4- The stability of the dynamic 



^A dynamic system x{t) = f(x{t — 1)) is additive hom.ogcneous of degree 1 if / is so, that is if Vx £ 
R", VA G M, f{X + x)^X + fix). 

^A dynamic system x{t) = f{x{t — 1)) is monotone if / is so, that is if Vxi, X2 £ K", aii < X2 => f{xi) < 
/(X2), where the order < is pointwise in M". 

■^A dynamic system x{t) = f{x{t — 1)) is non expansive if / is so, that is if there exists a norm || • || in R" 
such that Va;i,a;2 S R", ||/(a;2) - /(xi)|| < ||a;2-a;i||. 



system ([7]) is thus guaranteed from its non expansiveness. Moreover, ([7]) is connected (or 
communicating) □ jlll . Il2| . An important result from 18|, [ll|] (Theorem [1] below) permits 



the analysis of non expansive and connected dynamic systems. 



Theorem 1. llSAlAJ If a dynamic system x{t) = f{x{t — l)) is non expansive and connected, 
then the additive eigenvalue problem v + x = f{x) admits a solution {v, x), where x is defined 
up to an additive constant, not necessarily in a unique way, andv G M" is unique. Moreover, 
the dynamic system admits an average growth rate vector x, which is unique (independent 
of the initial condition) and given by xif) =^iv,v,- ■ ■ ,v). 

The model treated in this article can be seen as an extension of the min-plus model ([7]). 
In Section 2, we present the model. It is called piecewise linear car-following model because 
it is based on a piecewise linear fundamental diagram Ve- The model describes the traffic 
of cars on a ring road of one lane without passing. The stability conditions of the dynamic 
system describing the traffic are determined. Under those conditions, the car-dynamics are 
interpreted as a dynamic programming equation (DPE) associated to a stochastic optimal 
control problem of a Markov chain. The DPE is solved analytically. We show that the 
individual behavior law V^ supposed in the model is realized on the collective stationary 
regime. Finally, the effect of the stability condition on the shape of the fundamental diagram 
is shown. 

In section 3, we give equivalent results of those given in section 2 for the traffic on an 
"open" road (a highway stretch for example), and conclude with an example, where we 
simulate the transient traffic, basing on a piecewise linear approximation of the diagram ([6]). 



In Appendix B we give more details on the duality in traffic modeling, of using the 
functions n(t, x),x{n, t) and t(n, x), where t(n, x) denotes the time of passage of the nth. car 
by position x. 

2. Piecewise linear car following model 

The behavioral law Ve is an increasing curve bounded by the free speed vq. Moreover, 
Veiy) = for y G [0,yj] where yj denotes the jam inter-vehicular distance. We propose here 
to approximate the curve Ve with a piecewise-linear curve 

Ve{y) = mmmax{auwy + f^uw}, (10) 

ueu wew 

where U and W are two finite sets of indices. Since Vg is increasing, we have a„^ > 
0,V(u,w;) eUxW. 

We are interested here on the discrete-time first-order dynamics 

Xn(t + 1) = Xn(t) + min{a„(x„_i(t) - Xn(t)) + /3u}, (11) 



"'An additive homogeneous of degree 1 and monotone dynamic system x{t) = f{x{t — 1)) with x G M" is 
connected if its associated graph is strongly connected. The graph associated to that dynamic system is the 
graph with n nodes and whose arcs arc determined as foUows. There exists an arc from a node i to a node 
j if lim^^oofj{vei) ~ oo, where e^ denotes the i*'* vector of the canonic basis of K". 



and 



x„(t + 1) = Xnit) + minmax{a„^(x„_i(t) - x„(t)) + Puw}- 

u£U U)6>V 



(12) 



It is clear that (fT2l) extends ( fTTl) . The model ( fT2|) is also an extension of both linear model Q 
and min-plus model ([7]). 

In this article, we characterize the stability of the dynamics fll2p . calculate the stationary 
regimes, show that the fundamental diagrams are effectively realized at the stationary regime, 
and analyze the transient traffic. We will distinguish two cases: Traffic on a ring road and 
traffic on an "open" road. 

2.1. Traffic on a ring road 



We follow here the modeling of [27[. Let us consider v cars moving a one-lane ring road 
in one direction without passing. We assume that the cars have the same length that we 
take here as the unity of distance. The road is assumed to be of size n; that is, it can contain 
at most /i cars. The car density on the road is thus p = v/ p. 

Stochastic optimal control model 

We consider here the car dynamics (TTTj) . That is to say that each car n maximizes its 
velocity at time t under the constraints 



Xn{t + 1) < Xn{t) + au{Xn-l{t) - Xn{t)) + /?„, Vm G U . 



(13) 



Each constraint of ( [T3|l bounds the velocity Xn{t + 1) — x„(t) by a sum of a fixed term Pu 
and a term depending linearly on the inter-vehicular distance. 
Let us first notice that ( ITT]) , on the ring road, is written 



Xn{t + 1) = Xn{t) + min{a„(x„_i(t) - x„(t)) + /3„}, 



for n > 2, 



Xi(t) + min{a„(x„(t) + /U - Xi{t)) + /?„}, 



xi(t + i; 

which can also be written 

Xn{t + 1) = min{(l - au)xn{t) + a^Xn-iit) + 13^] 



for n > 2, 



xi{t + 1) = min{(l - a„)xi(t) + QuXnit) + a^i^/p + /3„}. 
Let us denote by M^, u alA the family of matrices defined by 



W 



\ — Olu 







a.u I- oiu 

and by c", u &U, the family of vectors defined by 



u t 



aui^/p + /3u, /3u, 



/3J. 



The dynamics ( ITT]) are then written : 

Xn{t+l) = mm{[M''x{t)]n + cl}, l<n<u. (14) 

ueu 

The system ( 1T^ is additive homogeneous of degree 1 by the definition of the matrices 
M", u eU. It is monotone under the condition that all the components of M", u (zU are 
non negative, which is equivalent to a„ G [0,1], Vn G U. Hence, under that condition, the 
system flT^ is non expansive. 

Moreover, the matrices M^,u G U are stochastic g. Those matrices can then be seen as 
transition matrices of a controlled Markov chain, where the set of controls is U. The connect- 
edness of the system flT^ . as defined in the previous section, is related to the irreducibility 
of the Markov chain with transition matrices M'^,u G U. It is easy to check that flT^ is 



connected if and only if 3u E lA^au ^ (0, 1]; see [Appendix A| for the proof. That condition 
is interpreted in term of traffic by saying that every car moves by taking into account the 
position of the car ahead. 

Consequently, under the condition Vn G W, a„ G [0,1] and 3u E U,au E (0,1], the 
dynamic system ^\M is non expansive and is connected. Therefore, by Theorem [H we 
conclude that the additive eigenvalue problem 

w + x„ = min{[M"x]„ + cJ^}, 1 < n < z/ , (15) 

describing the stationary regime of the dynamic system (TT4l) . admits a solution {v,x), where 
X is defined up to an additive constant, not necessarily in a unique way, and -u G M" is 
unique. Moreover, the dynamic system admits a unique average growth rate per time unit 
X, whose components are all equal and coincide with v. 

The average growth rate per time unit x of the system (fT4l) is interpreted in term of traffic 
as the stationary car-velocity. The additive eigenvector x gives the asymptotic distribution 
of cars on the ring, x is given up to an additive constant, since the car-dynamics (I14p is 
additive homogeneous of degree 1. That is to say that if {v,x) is a solution of (TT^ then 
(w, e + x) is also a solution for flTSl) . for every constant e G M. 

Let us now give an interpretation of the model in term of ergodic stochastic optimal 
control. Indeed, ( fT5l) can be seen as a dynamic programming equation of an ergodic stochastic 
optimal control problem of a Markov chain with transition matrices M",m G U and costs 
c**, u eU, and with a set of states A/" = {1, 2, ■ ■ ■ , u}. The stochastic optimal control problem 
of the chain is written 



I lim ^Y" C \ 



mmE{jim -y^c2}, (16) 

where 5 is a set of feedback strategies on A/". A strategy s E S associates to every state 
n E M a. control u ElA (that is Ut = s{nt)). 

The following result gives one solution [v, x) for the dynamic programming equation flTSj) . 



^We mean here M^- > 0, Vi, j and Y.j M^j = 1, Vi. 
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Theorem 2. The system l[TS\) admits a solution {v,x) given by: 

V =mm{auy + (3u}, 

x=*[0 y 2y ■■■ (u - l)y]. 



Proof. First, because of the symmetry of the system (1151) . it is natural that the asymptotic 
car-positions x„, 1 < n < z/ are uniformly distributed on the ring, and that the optimal 
strategy is independent of the state x. Let us prove it. 

Let it E U he defined by any + f3u = min„gi^{a„|/ + /?„} = v. Let x be the vector given 
in Theorem [2l Then Vn G {1, 2, ■ ■ ■ , z/} we have 



[M^x]n + cl = {any + /?«) + Xn = m.m{auy + /?«) + Xn = min{[M''x],, + c"} = t; + x„ 



D 



In term of traffic. Theorem [2] shows that the car-dynamics is stable under the condition 
«„ G [0,1], and the average car speed is given by the additive eigenvalue of the asymp- 
totic dynamics in the case where the system is connected. Moreover, it affirms that the 
fundamental diagram supposed in the model is realized at the stationary regime. 

V = mm{auy + /3n}, (17) 

q = min{a„ + /3„p}. (18) 

It is important to note here that, up to the assumption a„ G [0, 1], Vu G U, every concave 
curve Ve or Q^ can be approximated with fITTj) or flTSl) . Indeed, approximating fundamental 
diagrams using those formulas is nothing but computing Fenchel transforms; see (y, yj . More 
precisely, if we denote by V the set V = {/3u, u G U} and define the function g by: 

g: V ^ M 



then 



(1 = Qe{p) = mm{pv - g{v)) = g*{p), 



where g* denotes the Fenchel transform of g. 

Finally, we note that the min-plus linear model is a particular case of the model presented 
in this section, where U = {ui,U2} with {ai, (3i) = (0,w) and {a2, (32) = (1,— cr). In this 
case, the fundamental traffic diagram is approximated with a piecewise linear curve with two 
segments. 



Stochastic game model 

We consider in this section the car dynamics fll2p . again with the assumption \/{u,w) G 
U X yV,auw ^ [0,1] and 3{u,w) E U x yV,auw ^ (0,1]- The dynamic system (TT^ is 
interpreted as a stochastic dynamic programming equation associated to a stochastic game 
problem on a controlled Markov chain. As above, a generalized eigenvalue problem is solved. 
The extension we make here approximates non concave fundamental diagrams. 

In term of traffic, we take into account the drivers' behavior changing from low densities 
to high ones. The difference between these two situations is that in low densities, drivers, 
moving, or being able to m,ove with high velocities, they try to leave large safety distances 
between each other, so the safety distances are maximized; whilst in high densities, drivers, 
moving, or having to move with low velocities, they try to leave small safety distances 
between each other in order to avoid jams; so they minimize safety distances. 

To illustrate this idea, let us consider the following two dynamics of a given car n. 

Xnit + 1) = min{x„(t) + v, x„_i(t) - a}, (19) 

Xn{t + 1) = min{x„(t) + f , max{x„_i(t) - a, (x„(t) + a;„_i(t))/2}}. (20) 



The dynamics flT9|) is a min-plus dynamics which grossly tell that cars move with their 
desired velocity v at the fluid regime and they keep a safety distance a at the congested 
regime. The dynamics ( 12U]) distinguishes two situations at the congested regime: 

• In a relatively low density situation where the cars are separated by a distance that 
equals at least to 2o" we have 

max{x„_i(t) - (J, {xn{t) + x,„_i(t))/2} = x,„_i(t) - a. 

• In a high density situation, where the cars are separated by distances less than 2a we 
have 

max{x„,_i(t) - cr, (x„,(t) + x„_i(t))/2} = (x„,(t) + x„_i(t))/2. 

In this case, we accept the cars moving closer but by reducing the approach speed in 
order to avoid collisions. This is realistic. 

The situation we have considered in (1201) is realistic and very simple, but, it cannot be 
obtained without introducing a maximum operator in the dynamics (i.e. with only minimum 
operators). Indeed, with only minimum operators the approach is mechanically reduced 
with the increasing of the car-density (in fact this is the concaveness of the fundamental 
diagram). Because of the realness of such scenarios, we think that the fundamental traffic 
diagram should be composed of two parts, a concave part at the fluid regime, and a convex 
part at the congested regime. The dynamics flT2|) generalizes this idea. 

The dynamics flT2l) can be written 

x„(t + l) =minmax{[M""'x(t)]„ + c7}, l<n<u, (21) 



where 

^11.11) u * * * L« 



M" 






^ti.in J- ^1. 



and 

Similarly, we can easily check that the dynamic system ( 12T|) is non expansive under the 
condition auw ^ [0, 1], V(m, w) eUx W. Its stability is thus guaranteed under that condition. 
If, in addition, 3(m, w) eU x W, auw ^ (0, 1], then the system is connected. In this case, we 
get the same results as in Theorem HJ That is, the eigenvalue problem 

t; + x„ = minmax|[M"'"d„ + cri, l<n<u (22) 

admits a solution (■?;, x) given by: 

V = min max{auwy + /Suto}, 
x=*[0 y 2y ■■■ {u - l)y]. 

Moreover, the dynamic system fl2T|) admits a unique average growth rate vector x, whose 
components are all equal to v. 

The car-dynamics is then stable under the condition V(n, w) eU x W, auw £ [0, 1], and 
the average car speed is given by the additive eigenvalue of the asymptotic dynamics in the 
case where the system is connected (that is, if 3{u, w) eU x W, auw ^ (0, 1]). The behavior 
law supposed in the model is realized at the stationary regime. 

V = min max{auwy + Puw}, (23) 

q = minmax{a„u, + /3uwP}- (24) 

In term of stochastic optimal control, the system (l22l) can be seen as a dynamic program- 
ming equation associated to a stochastic game, with two players, on a Markov chain. The 
set of states of the chain is again M = {1,2,--- ,z/}. The chain is controlled by two players, 
a minimizer one with a finite set U of controls, and a maximizer one with a finite set W of 
controls. The transitions and the costs of the chain are given by the matrices M^^ and the 
vectors c"'", {u,w) eU x W defined above. 

The stochastic optimal control problem is 

minmax I^.^E | ^lini^ ^ J2 <« } ' (^5) 



where S is the set of strategies assoicating to every state n G A/" a couple of commands 
(m, w) E U X W. It is assumed here that the maximizer knows at each step the decision of 
the minimizer. 

We now give a consequence of the stabihty condition auw G [0, l],V(n,u;) G W x W, on 
the shape of the fundamental diagrams ( l23l) and ( 12^ . As shown in Figure [H, where we have 
drawn the fundamental diagram iQ (with vq = 14 meter by half second, and 7 = 7.5), the 
stability condition puts the curves fl23|) and f lM|) in specific respective regions in the plan. 
Indeed, for the diagram (l23l) . if we assume that Vg is bounded by Vq, K(i/) = 0, Vy G [0, i/j], 




20 40 60 80 100 

Inter-vehicular distance (meter) 




0.4 0.5 
car-density 



Figure 1: The effect of the stability condition auw G [0, 1] on the shape of the fundamental diagram. 

and that Vg is continuous (and increasing), then starting by the point (|/j,0), one cannot 
join any point above the line passing by {yj, 0) and having the slope 1, with any sequence of 
segments of slopes auw G [0, 1]. We can write 

Ve{y) < max(0, min(i;o, y - yo))- 

Similarly, on the diagram (^^, if we assume that Q^ is continuous and Qe{p) = 0, Vp G [pj, 1], 
then going back from the point (pj,0), one cannot attain any point above the line passing 
by {pj, 0) and (0, 1), with a sequence of segments having their ordinates at the origin (a„^) 
in [0, 1]. We can write 

Qe{p) < max(0, min(woP, 1 - p/pj))- 



3. Traffic on an open road 

We study in this section the traffic on an open road with one lane and without passing. 
We are interested in the following dynamics. 



xi(t + l)=xi(t)+t;i(t), 

x„(t + 1) = min„ge^ max„gw{x„(t) + a„^[x„_i(t) - x„(t)] + /3uw}- 



(26) 
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If we denote by M""* the matrices 



M"" 



1 ■■■ 

a„„, 1 - a„ 



and by c""'(t) the vectors 

for [u,w) eU X yV and t G N, then the dynamic system fj26l) is written 

Xn(t+1) =minmax{[M""'x(t)]„ + c7}, 1 < n < z/ . (27) 

It is easy to check that the dynamic system ( 127|) is additive homogeneous of degree 1, 
and is monotone under the condition \f{u,w) G W x W, a^^ G [0, 1]. Therefore, under that 
condition, fl27|) is non expansive. However, fl27|) is not connected for every (m, w) eU x W. 

We will be interested here, in particular, in the stationary regime, where vi{t) reaches a 
fixed value vi. For this case, the eigenvalue problem associated to f l26l) is 



f + Xi = Xi +Vi, 

v + Xn = min„6iv max^gwjxn + a„^[x,i_i - x„] + /3„„}. 



(2^ 



The system (1251) is also written 



w + x„ = minmax{[M""'x]„ + c7}, 1 < n < z/ , (29) 

where 

*-- [^1; Puwy ' ' ' ■> Puwi Puwj' 

Then we have the following result. 

Theorem 3. For all y E M. satisfying ramu^u laaxwewio^uwy + (3uw) = vi, the couple {v,x) 
is a solution for the system ([l^j, where v = Vi and x is given up to an additive constant by 

x = '[{n-l)y, {n-2)y, ■■■, y, 0]. (30) 

Proof. The proof is similar to that of Theorem[2J Let y G M satisfying min^^^^ inax^^y\![auwy+ 
Puw) = Vi. Let {u,w) eU X W such that au^^y + Puw = fi. Let x be given by (15U]) . Then 
Vri G {1, 2, ■ ■ ■ , I/} we have 

[M"'^x]„ + or = {au^y + {3^^,) + x^ 

= mmmax{auwy + /?««,) + x^ 

= minmax[M""'d„ + cr 

= Vi+ Xn- 

n 
11 



We can easily check that for {u,w) eU xW such that auw = and /3uw = vi, every inter- 
vehicular distance y E M satisfies the condition min„g2^max^gvv;(atj^y + /?„„,) = Vi. Thus, 
such couples [u, w) do not count for that condition. Theorem (IHl) can then be announced 
differently. Let us denote by Wu for u eU the family of index sets 

W„ = {iv e W, (a„», Am) ^ (0,t;i)}, 
and by y the asymptotic average car-inter-vehicular distance: 

. "^1 Puw /0 1 N 

y = max mm , [SI) 

u€U weWu auw 

where we use the convention a/0 = -|-oo if a > and a/0 = — oo if a < 0. Then Theorem [3] 
tells simply that if j^ G M (i.e. — oo < y < +C)o), then the dynamic system (!28|) admits a 
solution {v, x) where -u = t;i is unique, and where x is not necessarily unique and is given up 
to an additive constant by 

x = '[{n-l)y, {n-2)y, •••, y, 0]. 

Non-uniform asymptotic car-distributions can also be obtained. Let us clarify the follow- 
ing three cases. 

• If 3m G U,\/w G Wu, auw = and f3uw < vi, then y = -|-oo. In this case, the distance 
between the first car and the other cars increases over time and goes to -|-oo. The 
asymptotic car-distribution on the road is not uniform. 

• If Vm G W, 3w G W„, auw = and /3uw > vi, then y = — oo. In this case, the first car 
is passed by all other cars, and the distance between the first car and the other cars 
increases over time and goes to -|-oo. The asymptotic car-distribution on the road is 
not uniform. 

• If Vn G U^'iw G Wu,a„u, = and if vaSxiu&A'^^'y^w&Wu Puw = ^'i, then for all x G W, 
iyi.x) is a solution for the system (128|) . In this case, every distribution of the cars 
moving all with the constant velocity f i is stationary. 

The formula fl3T]) is the fundamental traffic diagram expressing the average inter- vehicular 

distance as a function of the car speed at the stationary regime. In the case where only a 

minimum operator is used in fl2B]) . the formula (1311) is reduced to the convex fundamental 

diagram 

vi - (3u .„„^ 

y = max . [62) 

ugu au 

Example 1 . In order to understand the transient traffic, let us simulate the car-dynamics fl2H]) . 
We take as the time unit half a second (1/2 s), and as the distance unit 1 meter (m). The 
parameters of the model are determined by approximating the behavior law (Q, with a free 



velocity vq = 14 m/ l/2s (which is about 100 km/h) and 7 = 7.5 as in 25|; see Figure |2l 
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Figure 2: Approximation of the behavioral law ^ with a piecewise linear curve. 



The behavior law is approximated by the following piecewise linear curve of six segments. 
V{y) = max{ai2/ + /3i, min{a22/ + 1^2, a^y + /^s, «4y + I^a, a^y + {3^, a^y + /^e}}, 
where the parameters at and /3j for z = 1, 2, ■ ■ ■ , 6 are given by 
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We simulate the piecewise linear car-following model associated to the approximation 
above. 

Xi{t) =Xi{t-l)+Vi{t), 

Xn{t) = Xnit - 1) + V{Xn-l{t - 1) - Xnit - 1)). 



(33) 



The velocity of the first car vi{t),t > is varied in the time interval [0, 1000], then fixed 
to the free velocity Vq = 14 m/ l/2s in the time interval [1000,3000], and finally fixed on 
a velocity that exceeds vq in the remaining time [3000,7200]. The average inter- vehicular 
distance is then computed at every time t, and the results are shown in Figure [3l 

The simple simulation we made here permits to have an idea of the traffic in the transient 
regime. Figure |3] shows how the average of the inter- vehicular distance is changed due to 
a changing in the velocity of the first car. In the right side of Figure [3l we compare the 
fundamental law assumed in the model with the diagram giving the average inter-vehicular 
distance (with respect to the number of cars) function of the velocity of the first car (a 
kind of macroscopic fundamental diagram). We observe that loops are obtained on that 
diagram in the transient traffic. The loops are interpreted by the fact that once the velocity 
of the first car is temporarily stationary, the velocities of the following cars, and thus also 
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Figure 3: Simulation results. On the left-side: the first car velocity (solid line), and the average inter- vehicular 
distance (dash line) functions of time. On the right side: the approximation of the behavior law ([5]) (dash 
line) , and the average inter- vehicular distance obtained by simulation in function of the velocity of the first 
car (solid line). 



the average velocity of the cars, make some time to attain the first car velocity. It can also 
be interpreted by saying that even though the cars have, individually, the same response 
to a changing in inter-vehicular distance; their collective response depends on whether the 
inter-vehicular distance is increasing or decreasing. The apparition of such loops is due to 
the reaction time of drivers. It can also be related to the number of anticipation cars in 
case of multi-anticipative modeling. However, one may measure on a given section, different 
car-flows for the same car density (or occupancy rate) depending on the traffic acceleration 



or deceleration, and interpret it as the hysteresis phenomenon [8|, |3l|, l32|, |15 



4. Conclusion 

We proposed in this article a car-following model which extends the linear car-following 
model as well as the min-plus model. The stability and the stationary regimes of the model 
are characterized thanks to a variational formulation of the car-dynamics. This formulation, 
although already made with continuous-time models, it permits to clarify the stimulus- 
response process in microscopic discrete-time traffic models, and to interpret it in term of 
stochastic optimal control. Among the important questions to be treated in the future, the 
impacts of heterogeneity and anticipation in driving, on the transient and stationary traffic 
regimes, based on the model proposed in this article. 

Appendix A. Connectedness of system ( 1141) 

Let X e W . We denote by /i : M'^ — )■ M'^ the operator defined by h{x{t)) = x{t+ 1), where 

That is 



Xn{t + 1), for 1 < n < u are given by the definition of the system ([T 

Xnit + l) = mm{[M''x{t)]n + cl}, l<n<i^ 
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• If 3m G W, such that a„ G (0, 1], then for all n G {1,2, ■ ■ ■ ,z/}, there exists an arc, 
on the graph associated to /i, going from n — 1 to n (n being cyclic in {1, 2, ■ ■ ■ , z/}). 
Indeed, 

Xn{t + 1) = (1 - au)Xn{t) + auXn-l{t) + /3n- 

Then since a„ > 0, we get: 

lim /i„(i/en-i) = lim [a„i/ + (5u] = oo. 

where Cn-i denotes the [n — ly^ vector of the canonic basis of W. Therefore the graph 
associated to h is strongly connected. 

• If Vm G U, a-a = 0, then we can easily check that all arcs of the graph associated to h 
are loops. Hence that graph is not strongly connected. 

Appendix B. Duality in traffic modeling 

We show in this appendix the duality in traffic modeling of using the three functions 

• n{t, x): cumulated number of cars passed through position x from time up to time t. 

• x{n,t): position of car n at time t (or cumulated traveled distance of car n from time 
up to time t). 

• t{n,x): the time that car n passes by position x. 

We base on the Lagrangian traffic descriptions given in the introduction, and give the equiv- 
alent traffic descriptions by using the functions n(t,x) and t{n,x), in both cases of discrete- 
time and continuous-time modeling. 

1. In Eulerian traffic descriptions, the function n{t,x) is used. The partial derivative 
dtn{t,x) expresses the car-flow q{t,x) at time t and position x, while —dxn{t,x) ex- 
presses the car-density p{t, x) at time t and position x. The equality dtxn{t, x) = 
dxtn{t, x) gives the car conservative law: 

dtk{t,x) + dxq{t,x) = 0. (B.l) 

The first order traffic model LWR [26|, l29| supposes the existence of a fundamental 
diagram of traffic giving the car-flow g as a function of the car-density p at the station- 
ary regime, through a function Qe, and that the diagram also holds for the transient 
traffic : 

g(t,x)=ge(p(t,x)). (B.2) 

Then fIB.ip and (IB. 21) give the well known LWR model 



dtp{t, x) + dxp{t, x)Q',{p) = 0. (B.3) 
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Also discrete-time and -space Eulerian traffic models exist. Those models are in general 
derived from Petri nets as in [2, llO, |9[. For example the traffic on a 1-lane road without 



passing can be described by 

n{t, x) = min{aa; + n{t - r^, x - 6x), a^+sx + n{t - f^+sx, x + 6x)}, (B.4) 

where 

a^ denotes the number of cars being in [x — 5x, x) at time zero. 

Tx denotes the free travel-time of a car from x — 6x to x. 

(ix denotes the space non occupied by cars in {x — 6x, x) at time zero. If we denote by 
Cx the maximum number of cars that can be in (x — 6x, x), then we have simply 

fx denotes the reaction time of drivers in (x — 6x,x). That is the time interval 
between the time when {x — 6x, x) is free of cars and the time when a car being 
in (x — 26x,x — 6x) starts moving to (x — 6x,x). If we denote by Tx the total 
traveling time (reaction time + moving time) of a car from (x — 26x, x — 6x) to 
{x — 6x,x), then we have simply fx = Tx — Tx- 

2. By using the variable t{n,x), the ffist order differentiation of t{n,x) with respect to 
n, denoted z is z{n,x) = —{t{n — l,x) — t(n, x)), while the derivative of t{n,x) in x, 
denoted r is r{n, x) = dxt{n, x). We notice here that z and r are interpreted respectively 
as the inverse flow and the inverse velocity of vehicles. A conservation law (of time) is 
then written 

dxz{n, x) + r{n, x) — r{n — 1, x) = 0. (B.5) 

The law (jB.Sp combined with the fundamental diagram r = Re{z) gives the model 

dxr{n,x) = R'{z{n,x))Ar{n,x), (B.6) 

where Ar(n,x) = r{n — l,x) — r(n,x). Note that, having a fundamental diagram 
V = Ve{q) giving the stationary velocity as a function of the stationary flow, the 
diagram R^ is nothing but Re{z) = l/Ve{l/z). 

Discrete-time-and-space modeling with the function t{n,x) also exist. The models are 
also inspired from Petri net, and dual dynamics to flB.4p are obtained. For example, 
using the same notations as in flB.4|) . the traffic on a 1-lane road without passing can 



be described by 

t{n,x) =meix{Tx + t{n-ax,x-6x), fx+sx + t{n + ax+sx,x + 6x)}. (B.7) 

Note here that a max operator is used rather than a min one. For more details on the 
duality of flB.4p and (IB.7P and the meanings in term of Petri nets, event graphs and 
min-plus or max-plus algebras, see |2|. 
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